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ABSTRACT 

We make use of a 500 ks Chandra HRC-S/LETG spectrum of the blazar 
H 2356-309, combined with a lower S/N (100 ks) pilot LETG spectrum of the 
same target, to search for the presence of warm-hot absorbing gas associated with 
two Large-Scale Structures (LSSs) crossed by this sightline, and to constrain its 
physical state and geometry. Strong (logN y// > 10 16 cm~ 2 ) OVII Ka absorption 
associated with a third LSS crossed by this line of sight (the Sculptor Wall, SW), 
at z = 0.03, has already been detected in a previous work. Here we focus on two 
additional prominent filamentary LSSs along the same line of sight, at z=0.062 
(the Pisces-Cetus Supercluster, PCS) and at z=0.128 (the "Farther Sculptor 
Wall", FSW). 

The combined LETG spectrum has a S/N of ~ 11.6 — 12.6 per resolution 
element in the 20 — 25 A, and an average 3a sensitivity to intervening OVII Ka 
absorption line equivalent widths of EW ovu ~ 14 mA in the available redshift 
range (z < 0.165). No statistically significant (i.e. > 3a) individual absorption 
is detected from any of the strong He- or H-like transitions of C, O and Ne (the 
most abundant metals in gas with solar-like composition) at the redshifts of the 
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PCS and FSW structures, and down to the above EW thresholds. However we 
are still able to constrain the physical and geometrical parameters of the putative 
absorbing gas associated with these structures, by performing joint spectral fit 
of various marginal detections and upper limits of the strongest expected lines 
with our self-consistent hybrid ionization WHIM spectral model. 

At the redshift of the PCS we identify a warm phase with logT = 5.35lo]i3 K 
and logNn = (19.1 ± 0.2) cm -2 possibly coexisting with a much hotter and sta- 
tistically less significant phase with logT = 6. K and logNn = 20. cm -2 
(l(j errors). These two separate physical phases are identified through, and 
mainly constrained by, CV Ka (warm phase) and OVIII Ka (hot phase) ab- 
sorption, with single line significances of 1.5 a each. 

For the second LSS, at z ~ 0.128, only one hot component is hinted in 
the data, through OVIII Ka (1.6c) and NelX Ka (1.2<r). For this system, we 
estimate logT = 6.6^ K and log Nh = 19.8±g;| cm" 2 . 

Our column density and temperature constraints on the warm-hot gaseous 
content of these two LSSs, combined with the measurements obtained for the hot 
gas permeating the SW, allow us to estimate the cumulative number density per 
unit redshifts of OVII WHIM absorbers at 3 different equivalent width thresholds 
of 0.4mA, 7mA and 25.8mA. This is consistent with expectations only at the 
very low end of EW thresholds, but exceed predictions at 7mA and 25.8 mA(by 
more than 2a). We also estimate the cosmological mass density of the WHIM 
based on the 4 absorbers we tentatively detect along this line of sight, obtaining 
qWhim _ (o.O21^QQ^g)(Z/Z ) _1 , consistent with the cosmological mass density 
of the intergalactic 'missing baryons' only if we assume high metallicities (Z ~ 
Z & ). 

Subject headings: intergalactic medium, large-scale structure of universe 



Introduction 



Cosmological hydrodynamical simulations predict the gradual formation of a local (z < 
1) web of low density (n^ = 10 -6 — 10 -5 cm -3 ), warm-hot (T = 10 5 — 10 7 K) intergalactic gas, 
connecting virialized halos (i.e. galaxies, galaxy groups and clusters of galaxies), permeating 
the large-scale structures (LSS) of which these systems are constituents , and ultimately 
providing the necessary fuel for their contin uous growth (ICen et al.l Il995t ICen fc Ostriker 
19991 ; iDave et al.l l2001t ICen fc Ostrikerl 120061 ) . This low-redshift intergalactic gas is largely 
the same primordial gas present at redshift higher than ~ 2 in a cool photo-ionized phase 
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(the so called 'Lyman-a Forest)', but in a much hotter and metal-enriched phase, because 
of efficient shock-heating during the continuous process of LSS assembly and growth in a 
non-linear Universe, and of strong feedback with the same structures for which it provides 
building blocks. Due to its high temperatures this IGM phase has been dubbed Warm-Hot 
Intergalactic Medium (WHIM). 

Electrons and baryons in the WHIM are shock-heated during their infall in the dark mat- 
ter LSS potential well, and settle in filamentary /sheet-like structures surrounding LSSs. Such 
matter is predicted to account for a sizable fraction (~ 50%) of all the baryons in the local 
(z < 1) universe, and it is thus considered the best candidate to host the baryons seen at high 
redshift and missing from the low redshift census (IFukugita et al.lll998t iFukugita &: Peebles 
20041 ). 



Given its high temperature the WHIM can only emit or absorb in the Far UV and soft 
X-rays, mainly through Li- through H-like metal transitions and bremsstrahlung continuum 
emission. However, at WHIM densities both line and continuum emission are highly de- 
pressed (due to the dependency of these mechanisms on the square power of the emitters 
volume densities), and on average well below the sensitivity of current instrumentations. 
Nonetheless, statistical techniques based on cross correlation of regions with excess diffuse X- 
ray emission in ROSAT, Chandra and XMM-Newton deep exposures, with large-scale galaxy 
distribution, have probably already allowed the marginal detection of the d e nsity peak of 



the W HIM distribution (e.g. IScharf et al.ll200d : IZappacosta et al.ll2002l . 120051 : Werner et al. 



20081 ). 



A far more promising way to detect the WHIM is through resonant absorption by highly 
ionized met als. Intervening WHIM filaments should imprint a 'forest', i.e. the so called X- 
ray Forest (IHellsten et al.l Il998l ; iPerna fc Loebl 119981 ) , by analogy with the HI Ly a Forest 
copiously seen at z > 2, of metal absorption lines onto the spectra of bright background 
sources, whose strength depends only linearly on the absorber density, and is therefore less 
suppressed than the corresponding emission. Predicted equivalent widths (EWs) from the 
most abundant ions in WHIM, range from 1 mA to < 20 mA for the densest filaments. The 
detection of even the strongest of such absorption lines (probing only the very high density 
tail of the WHIM distribution: overdensity 5 = n b / < n b >> 300, compared to the average 
density of the Universe: < n b >= 2 x 10~ 7 (1 + z) 3 (Q b h 2 / '0.02), where Q b is the baryonic 
density parameter and h is the Hubble constant in units of 100 km/s/Mpc) is extremely 
challenging with the limited sensitivity (A e ff < 60cm 2 ) and resolution (R=E/AE = 400 — 
800) of the current Chandra and XMM-Newton high resolution X-ray gratings. A > 3cr 
detection of an EW=20 mA absorption line, requires spectra of the background targets with 
S/N> 8 per 50 mA resolution element (i.e. > 70 net counts per bin). These can only be 
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obtained, for the brightest (~ 10 ergs 



cm 



soft X-ray (0.5 — 2keV) targets in the 



extragalactic sky (preferably blazars, because of their intrinsically featureless spectra), with 
> 0.5 Ms exposures, while in quiescence, or > 100 ks exposures in outburst. 

Such dense WHIM filaments are rare. According to hydro dynamical simulations, at 
z ~ one expects < 0.05 WHIM filaments with EW(OVII^ Q )> 20 mA, per unit redshift. 
The probability of havi ng one of such filaments along a random line of sight up to z = 0.3, 
is therefore only 1.5% (IGehrela I1986T ) . and only few z > 0.3 targets with quiescent F(0.5- 
2 keV)> 10~ n erg s _1 cm -2 are available (e.g. Conciatore et al., 2009, priv. comm.). 
The number density of OVII WHIM absorbers per unit redshift increases by almost two 
orders of magnitude by lowering the EW(OVII^ a ) detection threshold down to > 2 mA 
so dramatically increasing the probability of randomly detecting one of such filaments to 
P(z = 0.3)=97.5%. This observational strategy has the advantage to probe, in theory, 
the bulk of the WHIM mass distribution, but requires incredibly high S/N spectra (> 75 
per resolution element with the Chandra LETG, and > 180 per resolution element with 
the XMM- Newton RGS, for a > 3a detection), obtainable only with several Ms exposures 
on the brightest possible z > 0.3 targets in their quiescent states. Moreover, even when 
such high S/N spectra are obtained (for example by observing background target during 
exceptionally high outburst: iNicastro et al.l l2005al lbl) the clear assessment of the statistical 
significance of such weak 2-3 mA absorption lines is hampered by our limited knowledge of 
the in strument systematics, which is comparable t o the statistical unce rtainties on the lines 
EWs kaastra et alll2003l : Easmussen et al]l2007[ INicastro et al.l 120081 ) 



An alternative observational strategy is t o select lines of sight wh ere the probability to 



cross dense WHIM filament is enhanced (e.g. IZappacostal 12005 



200 6J). WHIM gas density 



and metallicity is predicted ( and possibly in t he UV observed IStocke et al.ll2006l ) to correlate 
with LSS galaxy overdensity ( jViel et al.ll2005l ). Thus, chances of intercepting a dense WHIM 
filament may be improved by carefully selecting bright sources in the background of extreme 
filamentary LSS co ncentrations ( Fig, [p . This observational strategy has been successfully 
exploited recently in lBuote et al. who detected a strong (EW~ 30mA) absorption line 

in two XMM- Newton RGS and Chandra LETG spectra of the blazar H 2356-309 (z=0.165), 
and tentatively identified it wit h OVII Kg at a redsh ift consistent with that spanned by the 
intervening Sculptor Wall (SW; Ida Costa et al.lll988l ). 



Here, and in our companion paper (IFang et al.l 120091 ). we report on the successful ex- 
tension of this observational program. H 2356-309 has been recently re-observed with the 
Chandra LETG for 500 ks, as part of an approved cycle 10 GO observational program. The 
main goal of this deep observation was to confirm, at higher significance, the SW OVII detec- 
tion (IFang et al.l 120091 ). Secondary objectives, were to confirm the presence of other (lower 



-5 - 



, „™ 0-5* 00h -0-50" , ™ 
1.00 s -p— -1.00 s 




Fig. 1. — Sky map and wedge diagram of the region of the Sculptor Wall where the blazar 
H 2356-309 is located. The upper sky map refer to the FSW and the lower to the PCS. 
Galaxies (black point) clusters and groups of galaxies (magenta circles) in the wall are taken 
from NED. The wedge diagram show galaxies inside the dashed blue box drawn in the sky 
map and report the sightline to the blazar as red dashed line. The galaxy and cluster catalogs 
contain objects belonging to different parent catalogs. Hence we point out that some holes 
in the projected galaxy distribution are artificially caused by this non homogeneity like the 
hole at north-east of the blazar position visible at RA = O.lhours and DEC = — 29.5deg in 
the PCS sky map. 
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significance) lines from the same absorber (Buote et al., 2010 in prep.). Here instead we focus 
on constraining the physical parameters of the putative WHIM gas content of two additional 
galaxy LSSs present along this line of sight, at z — 0.062 (the Pisces-Cetus Supercluster; 
PCS) and z = 0.128 (a farther wall which we will call Farther Sculptor Wall; FSW). In this 
paper we use all the existing Chandra LETG data of H 2356-309, to characterize the physical 
properties of the WHIM permeating these additional structures, and conclude by estimating 
the contribution of such dense gaseous component of LSSs, to the WHIM cosmological mass 
density. In SJ2] we describe the richness of LSSs along the line of sight to H 2356-309. In 
£J3] and §4] we present the data and describe their reduction and analysis. $5] is devoted to 
a critical discussion of our findings. In §6] we summarize our conclusions. Throughout the 
paper we adopt a A-CDM cosmology, with h = 0.71, Qm = 0.27, Qa = 0.73. 



2. The LSS Richness of the Line of Sight to H 2356-309 



Figure [Tj shows the wedge diagram of the line of sight to H 2356-309 in the declination 
range —33 < 5 < —29. Galaxies (black point), clusters and groups of galaxies (magenta 
circles) shown in the diagrams of Fig. [TJ are extracted from a number of different and 
non-hom ogeneous catalogs a nd galaxy surveys (including the 2dF Galax y Redshift Survey - 
2dFGRS. ICo~lless et alJliooS and 6dF Galaxy Redshift Survey - 6dFGRS. Ijones et dJliooi l. 
and are the result of a general query to the Nasa/Ipac Extragalactic Database (NEDjfl As 
such, these diagrams do not represent complete flux-limited sample of the actual galaxy 
distribution along this line of sight. 

Several strong LSS concentrations are clearly visible, and cross the line of sight to 
H 2356-309 at, at least, three different average redshifts: < Z\ >= 0.03 (the SW), < z 2 >= 
0.062 (the PCS), and < z 3 >= 0.128 (the FSW). Both the PCS and the FSW LSSs are 
significantly larger than the SW, and are delimited in the wedge diagram of Fig. [Tj by cyan 
dashed arcs. The two 2D sky map projections on the right of the wedge diagram of Fig. [TJ 
show the RA versus DEC extent of these two LSSs. 



The PCS (Burns fc Batuskilll984l ; iTullv et alJll992h is one of the richest nearby (z<0.1) 
superclusters. It is clearly visible in the SPSS an d 2dF redshift surveys as a remarkable 
filament of galaxies (IPorter fc Raychaudhuryll2005l ). The structure intercepted by the line 
of sight to H 2356-309 is a long filament of galaxies located on the plane of the sky at 
z=0. 06-0. 063, within < IMpc from the projected blazar position. 



1 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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The FSW is a conspicuous wall of galaxies, which originates from the Sculptor Super- 
cluster at z=0.11, and stretches out to redshift z ~ 0.135 crossing the blazar sightline at 
Az = 0.127- 0.129. 



Observations and data reduction 



The blazar H 2356-309 has been observed twice with the Chandra LETG, as part of 
two different observational programs. A first 100 ks LETG obse rvation was performed in 
October 2007 and has been already analyzed in lBuote et al.l (120091 ) . A second, deeper LETG 
observation was carried out over the September-December 2008 period, through ten different 
pointings with exposures ranging from 15 ks to 100 ks, for a total of 496. 4ks. The aim of this 
observation was to secure, with a conservative flux of 1.0 x 10 -11 ergs -1 cm -2 (0.5-2.0 keV) 
a 5cr significant detection of an abs orber with a colum n density of at least 9 x 10 15 cm~ 2 
(this is the 90% lower limit found by lBuote et al.ll2009[ ) by means of a long 0.5 Ms non-ToO 
observation. Table [TJ shows the log of the observations. 

We mention that H 2356-309 has also been observed with XMM for 130 ks and this 



observation has been analyzed in iBuote et al.l (120090 . We will not use this observation in 
our analysis since: 1) it adds only 70 counts in 0.06A (the XMM-RGS FWHM) which are 
~ 1/4 the net counts all the Chandra data provide in 0.05A; 2) does not reach the wave- 
lengths of the CV transition which, as we will see in §4.21 characterizes the most important 
intervening system; 3) makes the analysis in §4.21 very complicated due to the presence of 
several instrumental features very close to the lines we are investigating. 



3.1. Chandra reduction 

Each Chandra observation was reduced with the latest version of the Chandra Interactive 
Analysis of Observation software (CIAO v. 4.1.2, CALDB v. 4.1.2), following the standard 
processing procedures outlined in the HRC-S/LETG Grating analysis threadj, and applying 
a new (and still not standard) filtering procedure^ on the level 1 event files. This allowed 
us to greatly reduce the number of background counts, compared to the standard pipeline 
procedure, while losing only a negligible percentage of source counts, and so greatly increasing 
the S/N of the background subtracted source spectrum. 



2 http://asc. harvard.edu/ciao/guides/gspec_hrcsletg. html 
3 http://asc. harvard.edu/contrib/letg/GainFilter/ 
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Table 1: Log of the Chandra observations of H 2356-309. 



ObsID 


Instrument 


Exposure 


Date of observation 


8120 


HRC-S/LETG 


96.49 


2007 


-10- 


11 


10497 


HRC-S/LETG 


53.93 


2008 


-09- 


19 


10498 


HRC-S/LETG 


77.88 


2008 


-09- 


22 


10499 


HRC-S/LETG 


58.69 


2008 


-09- 


29 


10500 


HRC-S/LETG 


16.15 


2008 


-12- 


25 


10577 


HRC-S/LETG 


81.72 


2008 


-09- 


17 


10761 


HRC-S/LETG 


42.15 


2008 


-09- 


27 


10762 


HRC-S/LETG 


35.17 


2008 


-09- 


25 


10764 


HRC-S/LETG 


100.4 


2008 


-10- 


10 


10840 


HRC-S/LETG 


15.14 


2008 


-12- 


23 


10841 


HRC-S/LETG 


15.15 


2008 


-12- 


28 




Fig. 2. — Two portions of the spectral regions considered in §4.11 The left panel shows the 
Chandra spectrum around the rest frame OVII line region while the right panel shows the 
CV region. Each of the possibly identified lines (modelled in the spectrum with a Gaussian 
line profile) is labelled with colors indicating different absorbing systems. The fitted redshift 
of each line is showed in Table |3] 
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For each observation, we produced background light curves and inspected them visually 
to filter out periods of background flares. Only few observations were affected by short and 
mild flaring periods and their screening removed a total of just 30 ks, leaving a cleaned total 
562 ks exposure. 

To maximize the throughput, we used the Chandra LETG in combination with the 
Chandra High Resolution Camera for Spectroscopy (HRC-S) dispersion detector. This has 
virtually no spectral resolution, which prevents the separation of the different spectral orders 
dispersed by the LETG. We extracted 'all-order' negative and positive source and background 
spectra with the CIAO tool tg_extract. Due to impossibility of separating HRC-LETG spec- 
tral orders, the spectral modeling of these spectra can only be performed by pre-folding 
the fitting models with the sum of the convolution products of the redistribution matrices 
(RMFs) and ancillary responses files (ARFs) belonging to the first N positive and nega- 
tive orders, respectively. Here N is a number that depends on the intensity of the source, 
its spectral energy distribution, and the possible presence of strong spectral absorption or 
emission features, and has to be carefully chosen on an observation by observation base. 
For our LETG observations of H 2356-309, we verified., that in no case high-energy photons 
from orders higher than 6 contaminate the low energy spectrum by more than 3% (0.06% 
at the OVII wavelength). We decided to conservatively build our final "all-order" response 
matrix by adding up positive and negative orders up to N=10. For each order, we used the 
ftools task MARFRMF0 to 'multiply' the RMF by its corresponding ARF. This produced 
10 positive and 10 negative normalized RMFs, which we then co-added with the ftools task 
ADDRMF§. 

Finally, to maximize the S/N of our spectra, for each observation we co-added negative 
and positive order source and background spectra (add_grating_orders task) and response 
files (ADDRMF). 

Because of our co-adding procedure, particular care must be devoted to account for 
the wavelength calibration inaccuracy. An updated version of the degap polynomial coeffi- 
cients based on empirical wavelength corrections from multiple HRC-S/LETG observations 
of Capella has been release since CALDB 3.20. This update allow the correction of the non- 
linearities in the HRC-S dispersion relation improving the uncertainties across the detector 

4 http://cxc. harvard.edu/ciao/threads/hrcsletg_orders/ 
5 http://heasarc. gsfc.nasa.gov/lheasoft/ftools/caldb/marfrmf.html 
6 http://heasarc. gsfc.nasa.gov/lheasoft/ftools/caldb/addrmf. html 
7 http://cxc. harvard.edu/cal/Letg/Hrc_disp/degap. html 
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from 0.014A to O.OIOA (RMS deviation). This uncertainty is lower at shorter wavelengths 
and higher at longer wavelengths. Since we have coadded positive and negative orders we 
should consider a larger uncertainty. The propagation of the RMS deviations gives an un- 
certainty of 0.014A. However to be conservative we adopt 20mA of uncertainty since the 
updated wavelength corrections have been applied to the rest frame position of the strongest 
soft X-ray metal electronic transitions and the interpolation of the correction may not be 
strictly valid to the position of blue- and red-shifted lines. 



4. Spectral Analysis 

We performed all our spectral analysis with the fitting package Sherpa of CIAO (v. 
4.1.2). T he statistics w e used for our fits is the data weighted x 2 with the Gehrels variance 



function ( lGehrelslll986l ). the default in Sherpa. We first checked for variability of the broad 
band 10-50A continuum (flux and spectral shape) of H 2356-309, between the 11 Chandra 
observations. We grouped each spectrum at a minimum of 20 counts per bin, and modeled 
each data-set independently with a po wer law attenuated b y the sightline Galactic column 
of neutral gas (N H = 1.44 x 10 20 cm -2 ; Kalberla et al.ll2005l ). We found that both the source 



spectral shape and flux varied only moderately (12% and 20%, respectively) between the 
ten 2008 observations: the 0.3-1.0 keV flux ranges within the 1.25 — 1.5 x \§~ xl ergs~ x cm~ 2 
interval, while the power law best fitting spectral indices vary between T = 1.96 — 2.20. The 
0.3-1 keV source flux during the 2007 observation, is 50-80% lower than during the 2009 
observations, but the best fitting power law spectral index of T = 2.19 ± 0.06 is still in the 
range measured during the 2008 observations. Since both the 10-50 A flux and spectral 
shape of the target varied only moderately between the different observations, and because 
we are interested in the search of narrow spectral features, which are not affected by broad- 
band spectral-shape variability, we decided to co-add the HRC-S/LETG spectra of all the 
observations, to increase the final S/N per resolution element. The resulting spectrum has 
~ 290 net source Counts per 50 mA Resolution Element (CPRE) and ~ 120—160 background 
CPRE at 22 A , and so a S/N~ 11..6 - 12.6 at 22 A. This gives a theoretical > la sensitivity 
to absorption lines with EW> 12 — 13 mA, at 22 A. We grouped the data to half the nominal 
FWHM of LETG spectra with 25 mA per bin for all the subsequent analysis. 



4.1. Search for Intervening Absorption Lines 



Our main goal is to constrain the warm-hot gaseous content of the two filamentary 
galaxy structures here identified as PCS and FSW. At the expected WHIM temperatures 



- 11 - 



Table 2: Strongest Metal Transitions in Gas at T=10 5 — 10 7 K 



Ion 


Wavelength 


Oscillator Strength 


CV Ka 


40.2678 


0.648 


CVI Ka 


33.7360 


0.416 


OVII Ka 


21.6019 


0.696 


OVIII Ka 


18.9689 


0.416 


NelX Ka 


13.4471 


0.724 


NeX Ka 


12.1339 


0.416 



the most intense absorption lines are the Ka resonant transitions from He- and H-like C, O 
and Ne. Table [2] lists the rest frame wavelengths and oscillator strengths of these transitions. 

A crucial condition for the search for unresolved absorption lines in intrinsically fea- 
tureless spectra, is the accurate subtraction of the local continuum. A visual inspection 
of the broad band, 10-50 A, residuals of the co-added LETG spectrum of H 2356-309, af- 
ter subtracting the best-fitting absorbed power law, clearly showed the presence of 1-3 A 
broad systematics, particularly near or at the wavelength of the main instrumental edges 
of CI and OI. In general, the residuals deviated significantly from zero, in both directions, 
in several spectral regions. This makes the assessment of the actual significance of a pos- 
sible absorption line at such wavelengths, difficult. For each of the two redshifted systems 
(the PCS and the FSW), we therefore decided to isolate four 3-6 A broad spectral regions, 
around the absorption lines of interest (Table 1), at the following rest-frame intervals: 11-14 
A (NeX, NelX), 18.5-25.0 A (OVII, OVII), 31.5-34.5 A (CVI), 39-42 A (CV). We fitted each 
of these spectral intervals with continuum models including third-order polynomia. We then 
inspected visually the residuals to look for possible relatively broad 0.5-1 A deviations (due 
to instrumental effects), and, when needed, added to the best-fit polynomium 0.5-1 A broad 
emission or absorption Gaussian, to improve the modelling of the local continua. We iterated 
this procedure until we obtained a satisfactory description of all the local continua, and a 
visual inspection of the residuals revealed no further deviations. 

For each of the transitions listed in Table 2, we then added to the local best-fitting 
continuum a negative-only Gaussian (i.e. with flux allowed to be only negative in the fit), 
with wavelength allowed to vary within the redshift intervals of the PCS and the FSW and 
FWHM frozen to 10 mA (unresolved in the LETG), and re- fitted the data. For each line 
we estimate a single line significance in standard deviations, as the ratio between its best- 
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fitting flux, and its la error (computed with the Sherpa routine projection, by leaving both 
the continua and Gaussian normalizations free to vary). Four and three of the transitions 
listed in Table 2 were preliminarly, and tentatively, identified at a single line significance 
> la, for the PCS (CV, OVII, OVIII, and NelX) and the FSW (CVI, OVIII and NelX), 
respectively. For the remaining transitions we list their 3a EW upper limits. For comparison, 
and completeness, Table 3 also lists the result of our fitting procedure for the lines listed 
in Table 2, at z ~ 0. five of these transitions are detected at > la at z ~ (CV, OVII, 
OVIII, NelX and NeX). What discussed above provides only a preliminary attempt to check 
for the presence of the expected absorption lines at the redshifts of the Large-Scale galaxy 
Structures present along the lines of sight. The presence (or absence) of a given line (at the 
quoted significance level: Table [3]) is only used as a guidance for the detailed global fitting 
procedure presented in next section, and that that makes use of our self-consistent WHIM 
collisional ionization plus photoionization model. 

Figure [2] shows two spectral portions of the total Chandra LETG spectrum of H 2356- 
309, containing the three lines tentatively identified at the redshifts of the PCS and the FSW. 



4.2. Constraining the Physics of the Putative PCS and FSW Absorbers with 

a Self -Consistent WHIM Model 

In §4.11 we generally searched for the presence of the strongest absorption metal lines 
expected from highly ionized gas with temperatures in the broad interval T = 10 5 — 10 7 K. 
However, these transitions belong to ions with quite different ionization potential (e.g CV 
and NeX) and whose relative abundances critically depend on the ionization mechanisms at 
work, and on the actual gas temperature. For example, in collisionally ionized gas, OVII, 
being H-like and so quite stable, is virtually the only abundant ion of O within a broad 
interval of temperatures, within logT(K)=5.5-6.2. In the same temperature interval the 
lower ionization ions CV and CVI are relatively abundant, with CV decreasing monotonically 
from ~ 75% down to ~ 5%, and CVI raising from ~ 25% up to ~ 50% at its logT(K)=6 
peak temperature, and then decreasing again down to ~ 20%. On the contrary the higher 
ionization ions OVIII, NelX and NeX are only important at the high-temperature extreme 
of this interval, with relative abundances steeply raising from very low values up to 23% 
(OVIII), 96% (NelX) and 3% (NeX). Modeling the spectra of extragalactic sources crossing 
regions of the Universe with large galaxy overdensities (expected to trace WHIM filaments) 



http://cxc.harvard.edu/sherpa/ahelp/projection.py.html 
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with self-consistent ionized absorber models, may therefore provide useful constraints on the 
ionization state and column density of the putative absorbers embedded in these LSSs, even 
if the spectral signature of the gas are individually marginally detected (and/or upper limits 
are obtained) as long as they are modelled jointly 

Here we make use of an adaptation of the Photoio nized Absorber Spec tral Engine 



(PHASE; Krongold et al., 2003) code, for WHIM gas (e.g. iNicastro et all 120091 ). The code 
includes more than 3000 electronic resonant transitions (including metal inner-shell) from 
all elements lighter than Ni, and computes, for a given H equivalent column density, tem- 
perature and turbulence velocity of the absorber, the Voigt-profile folded opacity of each 
transition. The ionization balance in the gas is computed by perturbating the equilibrium 
collisional ionization balance at a given temperature T, with photoionization by the meta- 
galactic UV and X-ray background at a given redshift (the redshift of the absorber). Such 
second-order photoionization contribution depends uniquely on the baryon density n& in the 
gas (the lower the density the higher the contribution of photoionization), and it starts to 
be effective at n^l0~ 5 cm -3 . 

In our fitting procedure, we use the same initial methodology adopted in §4.11 For 
each of the two super-structures, the PCS and the FSW, we fit the 4 different narrow-band 
portions of the continuum where the main transitions lie, independently. For each spectral 
interval, our fitting model includes the best-fitting continuum determined in §4. II attenuated 
by our hybrid-ionization absorption WHIM model. For each spectral interval we leave free 
to vary in the fit, the continuum normalization and two out of the five parameters of the 
WHIM model, namely: the equivalent H column density N# and the temperature T of 
the gas. Both N# and T are linked to their same respective values in the 4 independently 
fitted spectral regions. The remaining parameters of our WHIM model are the turbulence 
velocity v (summed in quadrature to the thermal Doppler parameter of a given transition), 
the redshift z and the baryon density rib of the absorbers. The baryon density rib is highly 
degenerate with the electron temperature (which is set mainly by collisions in shock-heated 
WHIM gas), and modifies only slightly the ionization balance of the gas, compared with 
pure collisional equilibrium at a given temperature. Consequently, nb can be constrained 
independently of temperature only with data where several transitions from different ions of 
the same element are clearly detected. In our fit of the putative absorbers at the redshifts 
of the PCS and the FSW, we therefore freeze the gas baryon density to a typical WHIM 
value of n& = 10~ 5 cm -3 . Analogously, the turbulent velocity of the absorber (degenerate 
with the ion column density for saturated lines), can only be properly constrained when the 
single absorption lines are resolved and their profile clearly detected at high significance. We 
freeze this parameter to v — 100 km s _1 , comp arable to typical ya lues of Doppler velocities 



inferred by hydrodynamical simulations (e.g. ICen fc Fangl 120061 ) . Finally, for each of the 
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4 spectral regions, we first leave the redshift of the absorbers to vary independently over 
the entire redshift extent of the two super-structures (z = 0.060 — 0.063 for the PCS and 
z = 0.127 — 0.129 for the FSW), and then refine the fitting by freezing the redshift of the 
absorber in the spectral region where the most significant absorption line is detected to its 
best fitting value, and leaving the redshifts of the other absorbers in the three remaining 
spectral regions, free to vary between ±20 mA from the frozen redshift of the most significant 
absorption line, to account for the 90% systematic wavelength uncertainties in the HRC-S 
LETG. 

Table 4 summarizes the results of our fit. Errors are quoted at a 68% confidence limit. 
For each parameter listed in table 4, we compute errors by allowing the continuum normaliza- 
tion to vary within its la uncertainty and allowing the other free WHIM model parameters 
free to vary except in the case of the estimation of log Nh and log T errors where we fixed 
the redshifts to their best-fit value. 



4.2.1. The PCS filament 

For this structure, we found the possible co-existence of two distinct WHIM phases 
(Figure [3]): a warm phase, traced by CV absorption (Figure HJ), with logT = 5.35lo!i3 K 
and log Nh = (19.1 ± 0.2) x (Z/Z Q ) cm , and a much hotter and less constrained phase, 
traced by OVIII absorption (Figure E), with logT = 6.9±g"| K and log Nh = 20.1±^ x 
(Z/Z©) -1 cm- 2 . 

The redshift interval traced by the distribution of galaxies of the PCS, along the line of 
sight to H 2356-309, is 0.060 < z < 0.063. The two WHIM phases tentatively identified here 
have best-fitting redshifts consistent with each-other, and with the PCS redshift interval, 
namely: z warm = 0.0623 ± 0.0005 and z hot = 0.063 ± 0.001. 

We note that the physical parameters of the putative warm component of the PCS are 
constrained much better than those of the hot component (Figure [3]). This is because the 
opacity of absorbing gas to light metal transitions decreases with increasing temperatures. 
At the best-fitting temperature of the hot component of the PCS, only residual opacity 
from highly ionized O, Ne and Fe is present (Figure EJ), and the two free model parameters, 
temperature and N#, becomes highly correlated (Figure El red, dotted lines), and therefore 
only poorly constrained. At the temperatures of the warm PCS phase, instead, several strong 
transitions from a number of abundant ions can still produce enough opacity in the data, 
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Fig. 3. — 68%, 90% and 95% temperature and equivalent H column density confidence 
contours for the two putative WHIM phases (warm - solid cyan - and hot - dotted red) 
permeating the PCS. 

which makes T and N# virtually uncorrelated (Figure [3j cyan solid lines ). This can be 
also seen in Figure HJ where we show two portions of the LETG spectrum of H 2356-309 
with the strongest absorption lines from the warm component of the PCS, superimposed 
to the best-fitting continuum (red line) plus warm WHIM model component (blue line). 
The right panel of Fig. H] shows the strongest line of the best-fitting warm- component of 
the PCS, the CV Ka at a redshifted wavelength of A = 42.776 A. However, several other, 
moderately strong, lines are predicted by this model at shorter wavelengths, and are shown 
in the left panel of Fig. HJ namely the outer shell OVII Ka absorption line at A = 22.95A 
(superimposed on the 01 instrumental edge), and the two inner shell Ka transitions from 
OV at A = 23.74A (A = 22.35A rest frame) and OIV at A = 24.17A (A = 22.75A rest 
frame). The data are consistent with the presence of these lines, which tightly constrain the 
temperature and column density of the warm component of the PCS. 

4.2.2. The FSW 

For the FSW, the LETG data are consistent with the presence of one hot WHIM com- 
ponent (Figure [3]), traced mainly by Ka OVIII and NelX absorption (Figure [7]). The best- 
fitting WHIM parameter of this component are logT = 6.6^q' 2 K and logNn = 19.8to'g cm" 2 , 



9 We caution that the opening of the contours at low temperatures is mostly likely a numerical artifact, 
due to the hitting of the low-temperature end of the model grid, during the error-estimate procedure. 
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Fig. 4. — Portions of the LETG spectrum of H 2356-309 showing the strongest absorption 
lines from the warm component of the PCS. The best-fitting continuum model is shown 
in solid thick red, while absorption lines from the warm WHIM component are shown in 
solid thin blue. The CV Ka absorption line at a redshifted wavelength of A = 42.776 A is 
the strongest line of the best-fitting warm-component of the PCS (right panel). The other 
strong lines predicted by the model are the outer shell OVII Ka absorption line at A = 22.95A 
(superimposed on the 01 instrumental edge), and the two inner shell Ka transitions from OV 
at A = 23.74A (A = 22.35A rest frame) and OIV at A = 24.17A (A = 22.75A rest frame; left 
panel). Also shown in the left panel are the z ~ OVII Ka and the atomic and molecular 
01 Ka lines. 



and, as for the hot component of the PCS, they are only poorly constrained, due to the high 
best-fitting temperature of the gas. The tighter lower limit on the temperature, compared 
to that of the hot phase of the PCS, is due to the inconsistency of the data with any OVII 
Ka absorption stronger than EW> 5.8 mA (3cr limit; Fig. third panel). . 

The redshift interval traced by the distribution of galaxies of the PCS, along the line of 
sight to H 2356-309, is 0.127 < z < 0.129. The WHIM phase tentatively identified here has 



10 At the redshift of the putative FSW WHIM component, the OVII Ka falls at A = 24.25 A, a region 
of the detector free of instrumental features, unlike the PCS case, where the OVII Ka falls at the exact 
wavelength of the instrumental OI edge. 
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Fig. 5. — Portions of the LETG spectrum of H 2356-309 showing the strongest absorption 
lines from the hot component of the PCS. The best-fitting continuum model is shown in 
solid thick red, while absorption lines from the warm WHIM component are shown in solid 
thin magenta. The OVIII Ka absorption line at a redshifted wavelength of A = 20.17 A 
is the strongest line of the best-fitting hot-component of the PCS (right panel). The other 
predicted lines are the NelX Ka absorption line at A = 14.29 A and two strong L shell 
transitions from FeXVIII (A = 15.09A, A = 14.20A rest frame) and FeXVII (A = 15.96A, 
A = 15.02A rest frame; left panel). 



best-fitting redshifts consistent at lcr with the low redshift end of the FSW interval, namely: 
zfsw-whim = 0.126 ± 0.001. 

5. Discussion 

The line of sight to the blazar H 2356-309 is extremely rich in galaxy large-scale filamen- 
tary structures (Figured]). Other than the SW at least two distinct filaments of galaxies cross 
this line of sight, at average redshifts < z\ >= 0.0615 (the PCS structure) and < z 2 >= 0.128 
(the FSW structure). These very rich structures have line of sight extensions of D\ = 12 A 
Mpc and D2 = 8 Mpc, respectively, implying Hubble flow velocity ranges of Av? = 890 km 
s _1 and Az>2 = 580 km s _1 . 
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Fig. 6. — 68%, 90% and 95% temperature and equivalent H column density confidence 
contours for the putative WHIM gas permeating the FSW. 



The thermal broadening of C or O lines in gas with T< 10 7 K is b < 85 km s 1 (C) 
or b < 74 km s -1 (O), negligible c ompared to turbulenc e induced by peculiar motion of 



the structures. In simulations (e.g. ICen fc Ostrikerl 120061 ) WHIM intrinsic (i.e. excluding 



Hubble flow broadening) turbulence are observed to be of ~ 100 km/s. If the PCS and the 
FSW galaxy structures were homogeneously embedded by WHIM, therefore, the Hubble- 
flow braodening would be by far the dominant broadening mechanism, and the gas would 
imprint metal absorption lines with FWHM~ 890 km s _1 and FWHM~ 580 km s _1 , for the 
PCS and the FSW respectively, easily resolved by the HRC-S/LETG (FWHM=750 km s" 1 
at 20 A and FWHM=375 km s" 1 at 40 A). 

On the contrary, unresolved O absoprtion lines in the HRC-LETG must imply that the 
gas is homogeneously spread over a limited portion of these two galaxy structures, extending 
not more than 4 Mpc along the line of sight (corresponding to Hubble-flow broadening of 
0.021A at 20 A, i.e. la of the LETG Line Spread Function) 



5.1. On the OVII Bearing WHIM or Galaxy Halo Gas 



Our total Chandra LETG spectrum of H 2356-309 is sensitive to absorption line EW> 14 
mA at 22 A, at > 3cr confidence level. For unsaturated lines (6^200 km s -1 ) this cor- 
responds to OVII column densities Nqvii — 3.4 x 10 15 (1 + z)~ x cm -2 . At overdensi- 
ties 5 ~ 50(1 + z)~ 3 and temperature T~ 10 6 K, typical of the b ulk (~ 50%) of the 
WHIM density-temperature distribution (e.g. ICen fc Ostrikerl 120061 ). and assuming ho- 
mogeneity, these columns correspond to line of sight extensions of the filament of D > 
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Fig. 7. — Portions of the LETG spectrum of H 2356-309 where the strongest absorption 
lines from the putative WHIM component of the FSW lie. The OVIII Ka absorption line 
at a redshifted wavelength of A = 21.36 A is the strongest line of the best-fitting FSW 
WHIM model (central panel). The other moderately strong predicted line is the NelX Ka 
absorption line at A = 15.14 A (top panel). In the third panel we also show the region where 
the OVII Ka absorption line of the best-fitting FSW WHIM model falls: the data are clearly 
incompatible with the presence of such a line. 
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N b /n b = \\ Nav„/(AnfavTT)Vnp )(Z/Z m ) 
8.5 x 10~ 4 



^ 1A(Z/Z _i q ) 1 Mpc (where we assumed Ao = 
Grevesse fc Anders! Il989l and a relative fraction of OVII fovn — 0.9). The 



extensions of the PCS and FSW galaxy structures, if entirely permeated by WHIM gas rep- 
resentative of the bulk of its density-temperature distribution, should then guarantee the 
high-sig nificance detection of strong OVII absorbers. This has proven to be true only for 
the SW (iFang et ajj 120091 1 which is the most extended structure along the line of sight with 
D = 16.6 Mpc. 

We did not detect strong (i.e. Nov// > 3.4 x 10 15 (1 + z)~ l cmT 2 at > 3a), unresolved 
(i.e. b < 285 km s^ 1 at lcr) OVII absorption along either of the other two LSSs crossing 
the line of sight to H 2356-309, the PCS and the FSW. These non-detections tell us that 
any ionized gas at typical WHIM temperatures (T~ 10 6 K) embedded in the 4 Mpc cores 
(i.e. imprinting unresolved O lines: see §5) of the line of sight extent of the PCS and the 
FSW structures must have overdensities 5 < 18(Z/Z Q 10 ) _1 (1 + z)~ 4 (where the additional 
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(1+2;) 3 term comes from the (1 + z) 3 redshift dependency of < n b >)• For both the PCS 
{5 < UiZ/Zo^)- 1 at > 3a), and the FSW {5 < 11 (Z / '^o.io) -1 at > 3a) these overdensities 
lie on the lower end of the predicted overdensity interval for the bulk of the WHIM, ranging 
between 5 ~ 5 — 50. This is the opposite of what intuitively expected (though with a large 
scatter, e.g. Viel et al., 2005), that richer galaxy LSSs potentially harbor denser WHIM 
filaments. 

Alternatively the OVII bearing gas in these structures, could spread over the large extent 
of the galaxy super-structure along the line of sight, and so produce broader and resolved 
shallow (hence more difficult to detect) lines in the LETG (see §5) , or could have metallicity 
lower than 10% Solar, in either cases relaxing the above limits on the gas overdensities. 
However, any denser WHIM gas at the typical T ~ 10 6 K temperature and with line of sight 
turbulence velocity lower than a few hundreds km s _1 , should have been detected at high 
significance in the LETG spectrum of H 2356-309. 

An alternative interpretation of the OVII absorption, is that it could be associated with 
the hot extended halo of a single galaxy with small line of sight impact parameter. In this sce- 
nario, if r is the galaxy spherical halo radius, the impact parameter d must be smaller than r, 
and the line of sight can only cross a section of the halo D = 2rsin(a), where 0<a<7r/2is 
the angle between the direction of the galaxy-line of sight impact parameter d and the halo ra- 
dius r in the direction of the interception of the halo external boundary with the line of sight. 
By averaging over < a < rr/2, we get < D >= Ar/n. With this assumptions the baryon 
density of the galaxy halo is given by n b = N b / 'D<[[N OV ii I '(A fovii))/ '{^d/^)]( z I ' z &)~ 1 ■ 

The non-detection of OVII bearing gas up to the Novn > 3.4 x 10 15 (1 + z)' 1 cm -2 3cr 
limit, allows us to estimate stringent upper limits on the densities of putative galaxy halo 
gas with T ~ 10 6 K intercepting the line of sight to H 2356-309 at the redshift s of the PCS 



and t he FSW structures. By using an impact parameters 200 < d < 300 kpc (IStocke et al. 
2006h . we obtain n b ~ (3 - 5) x KT^Z/Zo.io) -1 cm" 3 for both the PCS and the FSW 



which is lower than the values derived for our galaxy halo or extended Local Group gas (e.g. 



Rasmussen et al.l 120031 ; IWilliams et al.l 120051 ) . These values would be raised at most by a 



factor of 2 by accounting for non constant gas density profile^]. 



11 We assumed a j3 — model ( Cavaliere fc Fusco-Femianol 1976 ) gas density distribution with (3 parameter 
ranging between 0.6 and 0.9 and assumed core radius ranging from 1 kpc to an unrealistically high 400 kpc 
value. 
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5.2. The Ionized Gas Content of the PCS and FSW Galaxy Structures 

As discussed in the previous section, we do not detect significant amount of OVII- 
bearing (i.e. typical of the bulk WHIM temperature-density distribution) gas in either 
the PCS or the FSW galaxy super-structures. However, at the redshift of these structure, 
we do marginally detect a number of (individually low-significance) metal absorption lines, 
from either low-ionization (OVI, OV, CV) or high-ionization (OVIII, NelX) ions. Such 
ions populate gas with temperatures in the low- or high-end of the WHIM temperature 
distribution, containing roughly 27% and 23% of the predicted WHIM mass, respectively. 
The absorption lines hinted in the LETG spectrum of H 2356-309 at the redshifts of the 
PCS and the FSW, are all unresolved, which implies line of sight extensions of the absorber 
D<A Mpc. 

Despite the low statistical significance of each of these individual absorption lines, we 
were still able to constrain the main physical parameters of the absorbing gas, namely its 
equivalent H column density and temperature, by modelling the broad-band LETG data 
with our hybridly ionized WHIM gas model. We identify two different absorbing WHIM 
phases at the redshift of the PCS, with logT = 5.351^3, log Nh = 19.1 ± O.2(Z/Z )- 1 , and 
logT = 6.9lo;g, logN H = 20.ll?;7(Z/Z Q )- 1 , for the warm and the hot phases respectively. 
For the FSW, instead, only one hot phase is tentatively detected, with logT = 6.6ta2 and 
logN H = 19.8^(Z/Z )- 1 . 

We can infer baryon densities lower limits for these systems, by assuming that the 
absorbers are embedded in their galaxy superstructures, and have dense cores extending < 4 
Mpc along the line of sight (we note that for the two phases in the PCS, assuming their 
are physically separated, the total extent can be close to the entire line of sight extension 
of the PCS galaxy filament). By conservatively assuming the —la Nh value as N b for the 
two PCS WHIM phases we get n b (Warm) = N b (Warm)/D > 6.4 x 10 -7 (Z/Zq)- 1 cm" 3 
(5 > 2.7(Z/Z Q )- 1 ) tmdn b (hot) = N b (hot)/D > 2.0xl(T 7 (Z/Z Q )- 1 (5 > O.9(Z/Z )" 1 ), while 
for the hot phase of the FSW we obtain n b = N b /D > 8.1 x IQ-^Z/Zq)- 1 (5 > 2.8(Z/Z Q )- 1 ), 
all consistent with predicted WHIM overdensities. 



5.3. Number density of OVII absorbers along the line of sight to H 2356-309 

From our best-fitting WHIM model temperatures for the three putative WHIM ab- 
sorbers at the redshifts of the PCS and the FSW, we can infer the relative fraction of OVII 
in each phase. These are fofn° S = 0-15, fovn° S = 9 x 10 ~ 4 and fovn = °- 028 - These ion 
fractions correspond to columns of ^ovn° S = L6 ±a6 x 1Ql5 cm ~ 2 , ^ovn CS = ^-§-11 x 1Ql3 
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cm" 2 and N^Ji 



1.5+?;| x 10 15 cm" 2 , or to unsaturated OVII EWs of 7±J mA, 0.42±°;* 2 



mA and 7^g X mA, respectively. 

Assuming the best-fitting temperature and N# values as face values, the lowest OVII 
equivalent width that we are indirectly probing, through our hybrid-ionization WHIM code, 
along the line of sight to H 2356-309. is therefore EW(OVII)>0.4 mA, corresponding to the 
hot phase of the PCS. Combining our three tentative WHIM detections, with the WHIM 
detection associated with the SW, leaves us with 4 distinct WHIM systems with EW(OVII)> 
0.4 mA, along a Az = 0.165 path length, or a number density of EW(OVII)> 0.4 mA 
filaments of dN(EW > 0A)/dz = 24.2^ 1 1 r' 2 ( allowing; for the large la uncertainties due to 
the small number statistics, i.e. 
simulation predictions (e.g 



Gehrels 



Cen fc Fang 



986). This is fully consistent with hydrodynamical 



2006|) 



We can, in principle, also compute the cumulative number density of OVII filaments 
per unit redshifts, with EW larger than a given threshold, for two additional E W thresh- 
olds: > 7 mA (3 absorbers) and > 25.8 mAA (1 absorber: Fang et al.l 12009 ). We get 
dN(EW > 7)/dz = 18.2+^ 7 ' 7 and dN(EW > 25.8)/dz = 6.1™. We plot such derived 
number density per unit redshif t, in Figure p supe rimposed to the theoretical cumulative 
dN(> EWthresh) / dz curve from Cen Sz FangJ ( 2006 ). The data-point at 25.8 mA, exceeds 
the predictions by more than 2a. This could be due to a line of sight selection bias, or may 
reflect the possibility that one or more of the detected absorption do not actually arise in 
tenuous WHIM filaments, but in much denser galaxy halos of one (or more) components of 
the LSSs with line of sight impact parameter of ~ 200 — 300kpc. It may finally be that 
the most marginally detected systems, the Hot-PCS and the FSW, given their large uncer- 
tainties in temperature and column density may just be explained as occasional statistical 
fluctuation. In this case, considering only the SW and Warm-PCS systems, the number 
density of EW(OVII)> 7 mA filaments would be dN(EW > 7)/dz = 12.1±J?;? which within 
the large uncertainties is consistent with the predictions (see Figure |S] dashed data point). 



5.4. Cosmological Mass Density of the WHIM from the H 2356-309 Line of 

Sight Absorbers 

Finally, we derive the Cosmological mass density of the WHIM as measured along the 
line of sight to H 2356-309, by assuming that all 4 different phases seen at the redshifts of 
the three SW, PCS and FSW galaxy superstructures are intervening WHIM filaments. By 
propagating in quadrature the large errors associated with the equivalent H column den- 
sity measurements of the absorbers, and those intrinsic with cosmic variance low-number 
statistics, we obtain: VL^ HIM = (O.O2lt[J;^)(Z/Z )- 1 . This is consistent with the cos- 
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Fig. 8. — Cum ulative numb e r den sity of OVII absorbers per unit redshifts compared to the 
predictions by ICen fc Fans (120061 ). The dashed line is an extrapolation of the theoretical 
curve. The dashed data point is obtained by excluding the two most marginally detected 
systems: the hot-PCS and FSW. 



mological mass density of intergalactic Missing Baryons (£i™ lssm 9 M ~ 0.016 ± 0.005, e.g. 



Fukugita fc Peeblesll2004l ) in the local Universe, but assumes solar metallicities. Metallicities 
of < 10% Solar, would give an Q^ HIM in excess by > la compared to Q™ lssin 9- IGM ^ a g am 
probably indicating an intrinsic line of sight selection bias, or simply reflecting the possi- 
bility that one or more of the detected absorptions are not imprinted by tenuous WHIM 
filaments. Again as in previous section we excluded from the analysis the most uncertain 
systems Hot-PCS and FSW and estimated an Q™ HIM = (O.OO26l^^)(Z/Z )" 1 which is 
consistent with the predicted missing baryons by assuming 10% solar metallicity. 



6. Conclusions 



We have analyzed a long integration (600 ksec) Chandra HRC-S/LETG spectrum of the 
blazar H 2356-309 to search and characterize the WHIM in two large-scale galaxy structures 
along this line of sight and, more specifically, the "Pisces-Cetus Supercluster" (PCS) at 
z=0.062 and the "Farther Sculptor Wall" (FSW) at z=0.128. These are structures more 
prominent and more distant with respect to the "Sculptor Wall", at z=0.03, where WHIM 
was identified in previous works through the detection of the OVII absorption line. 
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Although we do not detect significant individual absorption lines in the PCS nor in the 
FSW, the joint analysis of the marginally detected lines (as well as of stringent upper limits) 
through a self-consistent hybrid ionization spectral model, allow us to constrain the physics 
of the WHIM in the two farther superstructures. The main results are summarized in the 
following: 



At the redshift of the PCS we identify two distinct phases: a warm phase, with logT = 
5.35^>3 07 K and logN H = (19.1 ± 0.2) cm , and a much hotter less significant phase, 
with logT = 6.9+"- 1 K and logN H = 20. lt°{ 3 7 cm" 2 (lex errors). 



At the redshift of the FSW only one hot component is hinted in the data, with logT 
6.6+"- 1 K and logN H = 19.8±g;| cm" 2 . 



Under the assumption that the absorbers are embedded in their galaxy superstructures 
having baryonic column densities Nf, equal to the —la Nh value, and have dense cores 
extending <4 Mpc along the line of sight, we can infer conservative lower limits on 
the baryons densities in the two systems for the three absorbers. More specifically, for 
the two PCS WHIM phases we get 5 > 2.7(Z/Z )~ 1 and 5 > 0.9(Z/Z Q )-\ while for 
the hot FSW phase we obtain 6 > 2.8(Z/Z ) _1 , all consistent with predicted WHIM 
overdensities. 

By combining the constraints on t he OVH absorber s in the PCS and in the FSW, with 



the previous detection in the SW (IFang et al.ll2009l ). we derive the cumulative number 



density of OVII absorbers per unit redshift, as a function of the EW(OVII). While 
at low equivalent widths (EW(OVII)>0.4 mA) the absorbers number density is fully 
consistent with the predictions of hydrodynamical simulations, at EW(OVII)>10 mA 
the inferred absorbers number density exceeds significantly the theoretical predictions. 
The latter finding may result from a line of sight selection bias, or may reflect the 
possibility that one or more of the detected absorptions do not arise in WHIM filaments 
but in galaxy halos. We considered also the possibility that the two most marginally 
detected systems (Hot-PCS and FSW) may be just statistical fluctuations and not 
absorbers. In this case the number density of OVII absorbers per unit redshift would 
agree with the predictions within the large uncertainties. 

Finally, by combining the measurements in all absorbers we derive a cosmological 
mass density of the WHIM of fi^HiM = (o.O21+°;^)(Z/Z )- 1 , consistent with the 
cosmological mass density of intergalactic missing baryons in the local universe. Yet, 
if the WHIM metallicities are <10% Solar, then the resulting n^ mM is significantly in 
excess of the missing baryon density, again possibly indicating a line of sight selection 
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bias, or reflecting the possibility that one or more absorbers are not associated with 
WHIM. The exclusion of the two most marginally detected systems may bring in 
agreement the f2^ HIM value if we assume O.1Z . 
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Table 3: Best-Fitting Absorption Line Parameters and Ids for the PCS and the FSW 



Wavelength 


EW 


Significance 61 


Identification 


Redshift fc 


(in A) 


(in mA) 


(in a) 






Pisces-Cetus Supercluster (PCS) 


42.770 ± 0.015 


19.6 ± 13.1 


1.5 


CV Ka 


0.0621 ±0.0004 


35.76 - 35.86 


< 5.5 


NA 


CVI Ka 


0.060 - 0.063 


22.970 ±0.015 


6.2 ±5.5 


1.1 


OVII Ka 


0.0633 ± 0.0008 


20.160 ±0.015 


7.3 ±4.8 


1.5 


OVIII Ka 


0.0628 ±0.0007 


14.280 ±0.015 


5.0 ±3.3 


1.5 


NelX Ka 


0.062 ±0.001 


12.86 - 12.90 


< 3.2 


NA 


NeX Ka 


0.060 - 0.063 


Farther Sculptor Wall (FSW) 


45.38 - 45.46 


< 6.1 


NA 


CV Ka 


0.127-0.129 


38.042 ±0.015 


10.4 ±8.5 


1.2 


CVI Ka 


0.1278 ±0.0004 


24.32 - 24.39 


< 5.8 


NA 


OVII Ka 


0.126 - 0.129 


21.396 ±0.015 


7.7 ±4.9 


1.6 


OVIII Ka 


0.1279 ±0.0008 


15.161 ±0.015 


4.0 ±3.3 


1.2 


NelX Ka 


0.127 ±0.001 


12.46 - 12.47 


< 3.1 


NA 


NeX Ka 


0.127-0.129 


z ~ 


40.274 ±0.015 


38.5 ± 12.3 


3.1 


CV Ka 


0.0002 ±0.0005 


33.72 - 33.75 


< 7.1 


NA 


CVI Ka 


±0.0004 


21.561 ±0.015 


8.1 ±5.0. 


1.6 


OVII Ka 


-0.0019 ±0.0007 


18.952 ±0.015 


4.6 ±3.6 


1.3 


OVIII Ka 


-0.0009 ± 0.0008 


13.449 ±0.015 


7.6 ±3.3 


2.3 


NelX Ka 


±0.001 


12.131 ±0.015 


7.7 ±2.9 


2.7 


NeX Ka 


±0.001 



a Single Line Significance in Standard Deviations, evaluated as the ratio between the best 
fitting EW and its la error (see text for details). 

6 The error is derived from the systematic la wavelength uncertainty of ±15 mA due to the 
non-linearity of the HRC-S LETG dispersion relationship. 
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Table 4: Redshift and physical parameters resulting from the fits 



logT 


\ogN H 


Redshift 


(in K) 


(in (Z/Z & )~i cm" 2 ) 





PCS (z=0.06-0.063: Warm Phase) 



5.35+^ 19.1 ±0.2 0.0623 ±0.0005 

PCS (z=0.06-0.063: Hot Phase) 
6.9±g;^ 20. ll^ 0.063 ±0.001 

FSW (z=0.127-0.129) 

6.6l{ji 19.8ig;| 0.126 ±0.001 



